1 



Disappearance of quasi-fermions in the strongly coupled plasma 
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o ■ 

' We non-perturbatively study the fermion spectrum in the chiral symmetric phase focus- 

, ing on the effects of in-medium corrections for gauge boson. The fermion spectrum is derived 

- by solving the Schwinger-Dyson equation (SDE) with ladder approximation on the real time 

axis. It is shown that the peak of the fermion spectral function is broadened by in-medium 
effects for gauge boson compared with the peak obtained with the tree gauge boson propa- 
gator. The peak becomes much broader as the value of the gauge coupling increases. This 
' broadening is caused by multiple scatterings of fermions and gauge bosons included through 

the non-perturbative resummation done by the SDE. In particular, the Landau damping of 
' (— I ' \ gauge boson propagator plays an important role in the broadening. Our results show no 

^ ^1 clear peak in the strong coupling region, implying the disappearance of quasi-fermions in the 

I ' strongly coupled plasma. This indicates that quasi-particle picture may be no longer valid 

CIh' in the strongly coupled QGP. 

D ■ 

' §1. Introduction 

In the quark-gluon plasma (QGP) quarks and gluons, existing as elementary 
' degrees of freedom, have different spectra from those for the non-interacting quarks 

0^ ■ and gluons. In the weak coupling limit {g <^ 1), employing the Hard Thermal Loop 

(HTL) approximation, it is shown^-*''^^ that light quarks (such as u and d) and gluons 
have poles with masses proportional to the temperature T and the gauge coupling g 
but no widths. 



> 



in 



o , 

Q^ , Recent experimental studies carried out at RHIC, ^ on the other hand, sug- 

' gests that the QGP near the phase transition is a strongly interacting system. This 

result seems consistent with the findings of recent studies employing lattice QCD, 
^ I which show that the lowest charmonium state survives at T higher than the critical 

' ■ temperature (Tc)^^ 

For the thermal mass of quark near Tc, several analyses show Mq^sLvk/T ~ 1 based 
on one-loop calculation taking the gluon condensate into account, Brueckner-type 
many-body scheme and data from the heavy quark potential in the lattice QCD,^) 
Nambu-Jona-Lasinio model. In the work done in Refs. 8) and 9), an analysis 
employing the Schwinger-Dyson equation (SDE) is carried out to show that the 
thermal mass of fermion in the very strong coupling region saturates at Mqu^rk/^ ~ 
1, which is consistent with the result by a quenched lattice QCD with two-pole 
fitting. 

It is important to study not only the thermal mass but also the entire spectral 
function for clarifying the properties of quarks in medium. In Ref. 6), the imaginary 
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part of the quark self-energy is shown to be on the order of the thermal mass. 
Analyses in Refs. 8) and 9) show that peaks of the spectral function become broad 
in the strong coupling region. The SDE used in Ref. 8) is solved with the tree gauge 
boson propagator, which shows the broadening of the peaks with increasing coupling. 
This is understood from the fact that the probability of gauge boson emission and 
absorption from a fermion increases rapidly with the coupling. 

In the present work, we investigate the effects of medium modified gauge bo- 
son on the fermion spectral function. We employ the HTL corrected gauge boson 
propagator to include the characteristic properties of collective excitations of the 
gauge boson, i.e. quasi-particle poles and the Landau damping. We will pay special 
attention to the effect of Landau damping and show that it gives further broadening 
owing to multiple scatterings with fermions and gauge bosons, resulting no clear 
peak in the spectral function in the strong coupling region. 

The contents of this paper are as follows. In section [2l we introduce the SDE 
on the real time axis as an iterative equation for the fermion spectral function.^) 
In section [3l the fermion spectral function obtained in this way are shown focusing 
on the effects of medium corrected gauge boson. Section d] is devoted to a brief 
summary and discussions. In appendices, we show the phase structure in our frame 
work, and dependence on the difference of the method for solving the SDE and the 
gauge fixing condition. 



§2. The Schwinger-Dyson equation at Finite T 



The SDE for the fermion propagator S in the imaginary time formalism is ex- 
pressed as 

^ f d^k 

= g^T I j^;^ltiS{iuJrn,k)j,y'D'^''{iu;n - iuJm,p- k), (2-1) 



m,=— oo 



where U and D are the fermion self-energy and a gauge boson propagator, respec- 
tively, and uJn = (2n + 1)ttT is the Matsubara frequency for the fermion. The SDE 
was employed to describe the spontaneous chiral symmetry breaking in the strong 
coupling gauge theories. ^^^^^^^ When the value of the gauge coupling g in Eq. (|2-1|) is 
larger than the critical value gc, the chiral symmetry is spontaneously broken at low 
T. When g is smaller than g^ on the other hand, the chiral symmetry is not broken 
even at T = 0. We summarize the phase structure in appendix |X1 In the following, 
we consider only the chiral symmetric phase to investigate the fermion spectrum. 

We start with formulating the SDE on the real time axis following Ref. 9). 
The Matsubara summation in Eq. (j2Tp can be performed by using the spectral 
representation for fermion and gauge boson propagators as follows: 

ltiPf{z,k)j^p'^{z','p-k) , (2-2) 
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where / (n) is the Fermi-Dirac (Bose-Einstein) distribution function, and pj and p'^ 
are the spectral functions for the fermion and the gauge boson introduced as 



S{iuJm,p) = — j dz 



Pf{z,k) 



ILOm - Z 

Vf^'^iiLOn - iio^,p- k) = - I dz ^^^B (^';P-k) (2-3) 

Then we replace the Matsubara frequency as iun po + ie to obtain the retarded 
fermion self-energy as 

^BiPo,P) = -9' [4^ [dz [ dz'^-^ J (1 - fiz) + n(z')) 

J (27r)^ J J po + te- z- z' ^ 

X -f^pf{z,k)j,p''^'{z',p-k)\. (2-4) 
Taking the imaginary parts of both sides, we have 

lTa.ER{po,p) = T:g^ j -^—^ j ^dko (1 - f{ko) + n(po - A;o)) 

l^iPf{h^k)lvP^B{PQ-ko,p-k) , (2-5) 

where A is the ultraviolet cutoff introduced for regularization. Using the dispersion 
relation, the real part of the retarded fermion self-energy is expressed as 

Rei;,(po,p-) = lp Tdp^ ^'f^^^"'^"^ , (2-6) 
J~A Po - PO 

where P denotes the principal integral. The full retarded fermion propagator is 
expressed with the retarded fermion self-energy Uji(pQ,p) as 

iSR{po,P) = , , I. • (2-7) 

The spectral function for full fermion propagator is in turn expressed with the re- 
tarded propagator as 

Pf{Po,p) = -Im [iSR{po,p)] , (2-8) 
vr 

where p = \p\ and we remark that the fermion spectral function pf has spinor 
structure. Equations (j2-5p . (j2-6p . ()2-7p and (j2-8p make a closed form for the spectral 
function of fermion pf when the spectral function of gauge boson p'^ is a known 
function. In the method used in the previous work,^) we computed both real and 
imaginary parts of fermion propagator. The present formulation has the advantage 
of reducing the computation task, because we only compute the imaginary part of 
fermion self-energy from the loop integral, while the real part is computed from the 
imaginary part using the dispersion relation. 
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To make Eqs. (l2-5p - (l2-8l) a closed form, we choose an appropriate gauge boson 
propagator for the system. In the present work, to study the effects of collective 
excitations of gauge boson, we use the HTL corrected gauge boson propagator in 
the Coulomb gauge. In the previous work,*) we used the Feynman gauge to avoid 
double pole term. When we take account of the HTL correction to the gauge boson, 
we encounter the double pole term in the Landau gauge as well as in the Feynman 
gauge. In the present work, for simplicity of computation, we choose the Coulomb 
gauge in which there is no double pole term. We will discuss the dependence of the 
fermion spectrum on the gauge choice in appendix [B] in the case of tree gauge boson 
propagator. 

The spectral function of the HTL corrected gauge boson in the Coulomb gauge 
is expressed as 



P^"(?o, q) = PTiqo, q)Pr + PUqo, qW'^S''^, 



(2- 



where Pj^^ is the transverse projection = Pjl" = 0, and P^ = 6^^ — q^q-^ /q^- 
and p^^ are spectral functions for transverse and longitudinal parts of the HTL 
corrected gauge boson propagator. They are expressed as^^ 



PT{qo,q) = ZT{q)[5{qo - ajxiq)) - S{qo + u;T{q))] + PTiqo,q), 
PL{qo,q) = ZL{q)[S{qo - ^L{q)) - Kq^ + ^L{q))] +l3L{qo,q), 



(2-10) 
(2-11) 



where LOT{q) and ujiiq) are dispersion relations of the transverse and longitudinal 
gauge bosons, /3t and I3l are continuum parts of spectral functions and Zt and Zl 
are pole residues. They are expressed as 



Priqo, 



m'ix{l 



X 



2.6(1^2) 



1) 



Zriq) 



PL{qo,q) 
ujT{q){uJTiq) 



x^ + 



In 



x-\ 



m 



2x^(1 



+ 



g2 + 2m2 1 - f In 



x+l 


)] 


x-1 





+ TT^nigX^ 



Zdq) 



uJL{q){ujliq) 



Scoy^iq) - (coUq) - q^) ' "^"'^ q^{q^ + 2^2 - u^Hq)) ' 



(2-12) 



(2-13) 



(2-14) 



where rUg is the thermal mass of gauge boson and is related with the plasma frequency 
Up as iOp = ■\f2f?>mg = u;t,l(0). 

It is useful to parameterize the retarded fermion self-energy i7jj and the fermion 
spectral function as follows: 



Pf 



70 - r 



..p - 



-Im iSR{po,p) = ^ 



P 



P 



70 7/>+ + 7: 7oH 7/5 



P 



P 



(2-15) 



in which p+ (p_) is the fermion spectral function for (anti-) fermion sector. Sub- 
stituting Eqs. (|2-9p and (|2-15p into Eq. (j2-5p and making suitable projections, the 
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energy part (^^) and the vector part (i^]^.) of the retarded self-energy become 

ImU'}i{po,p) 



Im.S]i{po,p) 



where q^ = pQ - ko and q = \p- k\. Equations (12061) . OTTD . (IMj) . (I2^ and ([245]) 
make a set of iterative equations for the spectral functions p± in a closed form. 

§3. Numerical Results 

In the following, we fix the cutoff A as A/T = 5, and we show only p+ noting 
the relation p^{pQ,p) = p^{—pQ,p) obtained from the charge conjugation invariance. 

In Fig. [H we plot the fermion spectral function p-|_ for g = 1, 2 and 3, where 
three upper panels represent the spectral functions of fermion interacting with the 
tree gauge boson and three lower panels represent those of fermion interacting with 
the HTL corrected gauge boson with rUg = gTj^f^. In the upper panels, we see 
that /0+ has two peaks in the low momentum region, one in the positive energy 
region corresponding to the normal quasi- fermion, and another in the negative energy 
region corresponding to the anti-plasmino. Two peaks become broad as the coupling 
increases. In the higher momentum region, fermion spectrum approaches the one 
of free fermion. These results are consistent with those in the previous work.^^ 
Comparing fermion spectra in the lower panels with those in the upper panels, we 
find that two peaks in the lower panels are broader than those in the upper panels, i.e. 
two peaks become broader by the effects of the thermal correction to the gauge boson 
in medium. In lower panels, two peaks become broader as the coupling increases as 
in the upper panels. At 51 = 2, although we can see clear two peaks in the upper 
panel, two peaks in the lower panel are very broad, and the quasi-particle picture 
seems no longer valid. At 5 = 3 in lower panel there is no clear peak, which implies 
the disappearance of quasi-fermions. 

The broadening of the peak with increasing coupling, as we stated in the previous 
work,^) can be understood as the increasing probability of the gauge boson emission 
and absorption from a fermion. Because the SDE contains the effects of multiple 
scatterings with gauge bosons through the self-consistency condition, the peak ob- 
tained by the SDE is broader than the one obtained by the one-loop approximation 
(see Ref. 8)). 

The difference between the HTL corrected gauge boson and tree gauge boson 
lies in the difference of the value of mg/T: The HTL corrected gauge boson becomes 
tree gauge boson when nig/T = in Eqs. (|2T2p - (|2T4p . Here, we treat rUg/T 



■Kg 



(27r)3 

X [2pT(go,g) + PL(go,9)] , 
^ d^k '■^ 



^i. n f(i.\ , I \\ P+iko,k) + p-^iko,k) 
dko (1 - fiko) + n{qo)) 



ng 



(27r)3 
2pT{q,q)- 



dko (1 - /(/Co) + n{qo)) 
pk + (jp' + k'^) cos 9 — pk cos^ 



(2-16) 

P+{ko,k) - p-{ko,k) 



PLiqo,q) COS) 



(2-17) 
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Fig. 1. The fermion spectral function p+ in the Coulomb gauge. Fermion spectra in the upper 
panels are employing the tree gauge boson and those in the lower panels are employing the 
HTL corrected gauge boson with mg/T — g/ When we go from the left panels to the right 
panels, the coupling becomes stronger in both upper and lower panels: g = 1 (left panels), g = 2 
(middle panels) and <? = 3 (right panels). The figures are clipped at p+T — 2. 

as a parameter in Eqs. (|2-12p ~ (j2-14p and study mg/T dependence of the fermion 
spectrum. In Fig. [21 we plot the fermion spectral functions at p = for mg/T = 0, 
1/5, 1/3, g/V^ and 1 with g = 2 fixed. We see that two peaks become broad as 
the value of mg/T increases, and that two peaks merge to very broad one peak for 
mg/T = 1. (See the dot-dashed (cyan) curve in Fig. [2]) 

To understand the broadening of the peak with increasing value of mg/T, we 
consider the case in which the fermion interacts with HTL corrected gauge boson 
without continuum parts /3r and (3l in Eqs. (|2-10p and (j2Tip . We call this gauge 
boson "HTL pole" gauge boson: 

PTiqo,q)\u-^T^l^p^l^„ =ZTiq)[Siqo-UJT{q))-S{qo+ujTiq))] , (3-1) 
PLiqo,q)\^:-in^j^p„l^„ = ZL{q)[6iqo-ujLiq))-S{qo+ujL{q))] ■ (3-2) 

On the other hand, the HTL corrected gauge boson with continuum part is labeled 
as "full HTL" gauge boson: 

PT(go,g)|<.fuu HTL" = ^T{q)[S{qo - uJT{q)) - S{qo + uJT{q))] + PT{qo,q) , (3-3) 

Pliqo, 9)| "fuu HTL" 

= ZL{qmqo-iOL{q))-6iqo + iOLiq))]+(3Liqo,q) • (34) 

We plot the spectral functions of fermion at p = coupled with tree gauge boson, 
the "full HTL" gauge boson and the "HTL pole" gauge boson in Fig. [3j 
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-5 -4 -3 -2 -1 1 2 3 4 5 
Po/T 

Fig. 2. Fermion spectral functions for g = 2 at p — 0. The solid (red) curve is for = (tree), the 
dashed (green) for mg/T = 1/5, the short-dashed (blue) for mg/T = 1/3, the dotted (magenta) 
for mg/T = (?/\/6 and the dot-dashed (cyan) for mg/T — 1. 



0.7 p 
0.6 - 




Po/T 

Fig. 3. Fermion spectral functions for g = 2 and mg/T — 1/3 at p = 0. Solid (red) curve is the 
spectrum of fermion coupled with tree gauge boson, dashed (green) curve is coupled with "full 
HTL" gauge boson expressed in Eqs. (|3-3|) and p-4|) . and dotted (blue) curve is coupled with 
"HTL pole" gauge boson expressed in Eqs. (|3-ip and H3-2[) . 

The spectral function of fermion interacting with the "HTL pole" gauge boson 
shown by the dotted (blue) curve in Fig. [3]is sharper than that with tree gauge boson 
shown by the solid (red) curve. Because the SDE is a self-consistent equation, the 
internal fermion propagator has a thermal mass. The imaginary part of the fermion 
self-energy at one-loop, in which both the internal fermion and gauge boson have 
pole masses (denoted by mj and rUg), is forbidden in the energy region \mf — mg\ < 
\Pq\ < raf + rUg at p = owing to the energy-momentum conservation. While the 
imaginary part in the SDE is not completely forbidden but suppressed in the region, 
since the internal fermion in the SDE has a width. As was observed in a Yukawa 
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Fig. 4. A diagram which generates the imaginary part of the propagator for the fermion interacting 
with the HTL corrected gauge boson. 



model in Ref. 17), the suppression of the imaginary part leads to a sharp peak in 
the spectral function. As a result, the peak of the fermion spectrum interacting with 
the "HTL pole" gauge boson is sharper than that with the tree gauge boson. As 
nig increases, the above region becomes larger and the peak of fermion spectrum 
becomes even sharper. When the value of nig is very large, the fermion spectrum 
will approach that of the free fermion because the gauge boson decouples from the 
system. 

From the above analysis, we find that the continuum parts in Eqs. ()3-3p and 
(j3-4p play important rolls for the broadening of the fermion spectrum in the case of 
the "full HTL" gauge boson. Below, we give an interpretation for the broadening of 
the fermion spectrum with increasing nig/T: The dominant parts of Pt and (3l are 
proportional to as we see in Eqs. ()2-12p and ()2-13p . i.e. the Landau damping of 
gauge boson becomes larger as the value of nig/T increases. nig/T is proportional to 
the coupling at the point indicated by "A" in Fig. [H [Note that the coupling at "A" 
is different from g] Therefore the broadening can be interpreted as the increasing 
probability of scatterings of fermion and gauge boson at "A" . We would like to note 
that this effects is not included in quenched lattice simulations. 

At the last of this section, we consider the coupling dependence (Fig. [5]) and the 
cutoff dependence (Fig.[U]) of the peak position of the fermion spectral function which 
we call the thermal mass. Figure [5] shows that the thermal masses hardly depend 
on the value of nig/T when the peaks exist as clear ones. On the other hand, the 
widths strongly depend on nig/T as seen in Fig. [21 and there are no quasi- fermions 
in the strong coupling region as seen in Fig. [TJ That is why the curves for nig/T > 
in Fig. [5] terminate at some values of the coupling. This is contrasted to the previous 
work,*) in which tree gauge boson {nig/T = 0) was used. There, we saw that the 
peak position saturates at a value of coupling and the peak position M becomes 
M/T ~ 1 in the very strong coupling region. From Fig. El we see that M/T depends 
little on A/T in the region A/T > 4. We have also checked that the shape of the 
fermion spectral function around the peak hardly depends on the cutoff. 

§4. A summary and discussions 

In this paper, we investigated the fermion spectra in the chiral symmetric phase, 
focusing on the effect of medium correction for gauge boson, by employing the 
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0.5 1 1.5 2 2.5 3 3.5 4 

g 

Fig. 5. The coupling dependence of the thermal mass M/T. The solid (red) curve represents the 
one for rUg — (tree), the dashed (green) curve for nig/T = 1/3, the short-dashed (blue) curve 
for vfiglT = 1/2, the dotted (magenta) curve for nig/T = l/\/6 and the dashed (black) line 
represents the thermal mass of fermion in the HTL approximation. In the case of vrig/T = 1, 
we cannot read the peak as an excitation. 

0.8 I 1 1 1 1 1 1 1 1 1 
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0.3 - 



2 I I I I I I I I I I 

2 4 6 8 10 12 14 16 18 20 

A/T 

Fig. 6. Cutoff dependence of the thermal mass M/T. The solid (red) curve is for mg/T — 1/3 and 
the dashed (green) for mg/T — g/\f%. For A/T > 4, the ratio M/T is almost independent of 
the cutoff. 

Schwinger-Dyson equation (SDE) for fermion with the ladder approximation. 

We found that in-medium effects for gauge boson propagator make the peak of 
the fermion spectral function broader compared with the previous result*^ obtained 
with the tree gauge boson propagator. The broadening occurs owing to multiple 
fermion-fermion scatterings in addition to fermion-gauge boson scatterings, included 
through the Landau damping of in-medium corrected gauge boson. Thus, it is im- 
portant to take into account of not only the thermal mass but also the Landau 
damping of gauge boson propagator. Our results show no clear peak in the strong 
coupling region, which implies the disappearance of quasi-fermions in the strongly 
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coupled plasma. 

We would like to make a comment on the Van Hove singularities. They exist 
in the case of the self-energy consisting of two very sharp particles whose dispersion 
relations are different from each other as in the fermion and gauge boson in the 
HTL approximation. However, if the particles have widths, these singularities are 
weakened. In the SDE, because the (internal) fermion has a broad width by non- 
perturbative effects through the self-consistency condition, we cannot see the effects 
of Van Hove singularities in the present analysis. 

In this paper, we solved the SDE in the Minkowski space as an iterative equation 
for the fermion spectral function. In the previous work,^) on the other hand, the 
solution is obtained by performing the analytic continuation from the imaginary time 
axis to the real time axis through an integral equation. There is little difference in 
the fermion spectra obtained by two methods as seen in appendix [Bj We also show 
that there is little gauge dependence of the peak position in appendix [Bl 

We discuss the application of the results of the present work to QCD. In the 
very high temperature region of QCD, the HTL approximation is valid for studying 
the quark spectrum. In the HTL approximation, the fixed coupling, which is the 
running coupling at an energy scale of order T, say g{T), is used, because the quarks 
and gluons with energy T give the dominant contribution. In studying the quark 
spectrum in the lower T region, the diagrams other than the HTL start to give non- 
negligible corrections, since the (fixed) coupling g{T) is larger. For this reason, we 
can state that, in our approach, we include a part of non-perturbative corrections 
from a certain diagrams by solving the Schwinger-Dyson equation. However, the 
following points of QCD are not included: (1) the asymptotic freedom through the 
running effects of the gauge coupling, (2) vertex corrections and (3) higher oder 
corrections for the gauge boson. 

First, we discuss the effect of asymptotic freedom. In our analysis, the fixed 
coupling is used, which can be understood as the running coupling at an energy scale 
of oder T, i.e. g{T). In the energy scale E lower than T, the gauge coupling does not 
run, and we have g{E) ~ g{T). This is because quarks and gluons decouple in this 
energy region owing to the thermal mass. In ii^ > T, gauge coupling is running, and 
we have g{E) < g{T). The difference between g{E) and g{T) might cause different 
results when we use the running coupling in the SDE analysis. However, in the 
present analysis, we have seen that peak position only depend on T which is an 
infrared scale, but does not depend on the cutoff A which is an ultraviolet scale: The 
peak position is determined by the dynamics of the infrared scale only. Therefore, 
we think that the peak positions at least will hardly change from those in the present 
analysis. 

Second, we discuss vertex corrections. When we take account of some vertex 
corrections, the peak of fermion spectrum will become even broader than that in 
the present analysis owing to multiple scatterings. Then the following qualitative 
structure of the quark spectrum will hardly change: The peak position of the quark 
is proportional to the coupling for small g{T) and the quasi-particle picture will be 
even worse. 

Third, we discuss higher order corrections for the gauge boson propagator. There 
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is a coupled SDE for fermion and gauge boson, in which we can include certain higher 
order corrections for the gauge boson propagator. The peak of fermion spectrum will 
become even broader than that in the present analysis owing to multiple scatterings. 
The peak position does not depend on a detailed structure of gauge boson propagator 
in the present analysis, which indicates that the peak position does not change from 
the present value in the small coupling region, and that the quasi-particle picture 
will not be good in the strong coupling region. 

To summarize, the following qualitative structure of the quark spectrum will 
hold in hot QCD: In the high T region where g{T) is small, the thermal masses of 
the quasi-quark and the plasmino are proportional to the coupling. In the low T 
region where g{T) is large, the width will become broad due to multiple scatterings. 
Especially, in the strongly coupled QGP near Tc, the quasi-particle picture may be 
no longer valid. 

It is suggested that bound states for light quarks exist near Tc in Ref. 18) and 
19). It is also suggested that new mesonic states consisting of a quark and a plas- 
mino exist in Ref. 20). It is very interesting to study these bound states for light 
quarks through the Bethe-Salpeter equation using the fermion spectra obtained in 
the present analysis. 
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Appendix A 

The SDE in the imaginary time formalism and the chiral phase transition 



In this appendix, we determine the critical temperature Tc of the chiral phase 
transition from the ladder SDE. We solve the SDE in the imaginary time formalism 
because the numerical cost is much smaller than solving the SDE on the real time 
axis. 

The ladder SDE in the imaginary time formalism is given by Eq. (|2-ip . The full 
fermion propagator in the imaginary time formalism is restricted by the rotational 
invariance and the parity invariance as follows: 

S{iujn,p) = r—-rr- — TTT- V • (^■^) 

B{iLOn,P) + A{iiL>n,P)P ■ 7 - C(za;„,p)za;„7o 

The HTL corrected gauge boson propagator in the Coulomb gauge is expressed as 

V^''{iujn - iujm,p- k) = ^riii^n - i^miP- k)PJ/:'' + Aiiiu^n - iujm,p- /c)(5^°5''°, 

(A-2) 

where At and A l are transverse and longitudinal HTL corrected gauge boson prop- 
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agators. They are expressed as 



AtHuj, q) = ^, (A-3) 

Z\i(iw, q) = — ^. (A-4) 

Q2 + 2m2 l-Mf^_2tan-iM 



2q 

Note that Eq. (|A-3p and Eq. (|A-4p have no poles for real uj. Using these fermion 
propagator and gauge boson propagators, the ladder SDE (j2-ip is decomposed into 
three coupled equations for B, A and C as 



■m = — rv-i " V / 



(A-5) 



m=—oo 



^ = 1 + ^ V / 

m=—oo " ^ ' ^ ^ ^ 

(A.6) 



po ^ J (27r)3 52 + C2a;^ + A2fe2 

Because the integrals in Eqs. ()A-5p - ()A-7p are divergent, we introduce the three- 
dimensional ultraviolet cutoff A. We also truncate the infinite sum of the Matsubara 
frequency at a finite number A^. However, A^ can be taken to be sufficiently large so 
that the following results do not depend on it. The angle integral in Eqs. ()A-5p - ()A-7p 
are numerically performed. 

In Fig. [71 we plot the coupling dependence of Tc/ A for nig/ A = 0, 0.1 and 0.2. 
We see that Tc/ A decreases as rUg/ A increases owing to the screening effects. This 
is consistent with the results in Ref. 14) which show that becomes small by the 
effects of the Debye screening. 

Appendix B 

Dependences on the gauge fixing condition and method for solving the SDE 



There are two methods to obtain the solutions of the SDE on the real time 
axis. One is to perform the analytic continuation for the solutions of the SDE in 
the imaginary time formalism through an integral equation, 2^) which was used in 
the previous work.®) Another method is to solve a self-consistent equation as an 
integral equation for fermion spectral function.^) In this appendix, we show that the 
dependence of the peak position of the fermion spectral functions on two methods 
together with the gauge dependence of the peak position. In Fig. [HI we plot the peak 
positions of fermion spectral functions obtained by the method 1 in the Feynman 
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Fig. 7. Coupling dependence of the critical temperature Tc/A. The chiral symmetric (broken) phase 
is above (below) the line. The solid (red) line shows Tc/A for m/A = 0, the dashed (green) line 
for 7n/A = 0.1 and the dotted (blue) line for nig/ A — 0.2. 

gauge, by the method 2 in the Feynman gauge and by the method 2 in the Coulomb 
gauge. We find that the peak positions httle depend on neither the method nor 
the gauge fixing. This is consistent with the result of Ref. 22) in which the exact 
calculation of the fermion self-energy at one-loop level is performed and it is shown 
that the peak positions little depend on the gauge choice. 
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Fig. 8. Coupling dependence of the thermal mass. The solid (red) curve represent the thermal mass 
obtained by the method 1 in the Feynman gauge. The dashed (green) curve and dotted (blue) 
curve represent the thermal mass obtained by the method 2, in the Feynman and Coulomb gauge, 
respectively. The dashed (black) line represents the thermal mass in the HTL approximation. 
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